Growth of Oscillations of a Ray about 

the Irregularly Wavy Axis of 

a Lens Light Guide 

By D. W. BERREMAN 

(Manuscript received July 20, 1905) 

// a ray is launched in a direction coincident with the axis of a lens light 
guide whose axis is somewhat wary, the ray will soon begin to oscillate about 
the axis. The amplitude of these oscillations will grow in proportion to the 
square root of the product of the distance from the origin and the natural 
wave number of the oscillations, on the average. The growth rate is propor- 
tional to the amplitude of the components of the spectrum of the waves in 
the guide axis in the immediate vicinity of the natural wave number of the 
ray oscillations. The rest of the spectrum of waves in the guide does not 
contribute appreciably to the growth of ray oscillations after the first few 
oscillations. A tractable analytic expression with four adjustable parameters 
is used to approximate the wave spectrum of the shape of the guide axis. The 
expression is used to illustrate the relations between various factors such as 
mechanical stiffness of the guide, spectrum of external forces, over-all am- 
plitude of waves in the guide, etc., and the rate of growth of the ray oscilla- 
tions. Estimates of the order of magnitude of the oscillation growth rate in 
some realistic models of light guides are made from these relations. 

I. INTRODUCTION 

The possibility of guiding a beam of light over a long distance through 
a series of lenses or a continuous lens-like medium has recently received 
considerable attention because such a system might be useful in com- 
munication. One problem that has only recently come under investiga- 
tion is the statistical growth rate of the oscillation of an initially paraxial 
ray about the axis of such a guide when the axis is crooked in a some- 
what random but partially coherent way. 1,2 In this paper I will show 
that the ray will oscillate about the axis with ever-increasing amplitude, 
on the average. 

The oscillations of the ray about the axis are analogous to the os- 
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(•illations of an undamped mechanical oscillator subject to a particular 
steady "noise spectrum" of forces. Spherical aberration in the guide 
would be analogous to anharmonicity in the mechanical oscillator. 
Distance in the optical problem corresponds to time in the mechanical 
problem. 

The dominant term in the expression for the amplitude of ray oscilla- 
tions is proportional to the square root of the distance from the origin, 
if the amount and spectral form of the crookedness is constant and if 
the lenses have no spherical aberration. The oscillations of the ray about 
the axis have a natural wave number, k c . The growth rate of ray os- 
cillations is also proportional to the square root of k„ . The third factor 
governing the growth rate is the amount of crookedness of the axis. A 
major portion of this paper is devoted to a mathematical description of 
the amount of crookedness of the axis of the guide. Only those compo- 
nents of the crookedness having approximately the same wave number 
as the natural oscillations of the beam have an appreciable influence on 
the growth of the oscillations beyond the first few oscillations. 

Some numerical examples of growth rates of oscillations in guides 
having reasonable amounts of crookedness or waviness are given at the 
end of the paper and in Table I. The amount of crookedness in the 
examples was obtained by estimating the variation of forces on a pipe of 
reasonable stiffness lying on a rough surface. See Fig. 1. The pipe was 
assumed to be relatively straight before the irregular forces due to the 
rough surface were applied. This specific method of estimating crooked- 
ness does not limit the generality of the relations between the crooked- 
ness spectrum and the growth rate of the oscillations. 

The results indicate that, unless such a light guide could be kept ex- 
tremely straight or free of waves on a scale approximating the wave- 
length of natural oscillations of the beam, the beam would have to be 
recentered at frequent intervals* or the guide would have to have a very 
large aperture to avoid vignetting of an initially paraxial ray. 

II. EQUATIONS OF THE RAY TRAJECTORY 

In order to keep the analysis simple, I will suppose that the axis of 
the guide is wavy only in one dimension, y. The distance from a straight 

* There is no physical reason why we could not redirect the beam down the axis 
of the guide at intervals with almost no loss of beam energy. For example, the 
position and direction of the main beam could be located at one point by reflect- 
ing a very small fraction of it into a photoelectric analyzing device with a very 
weakly reflecting mirror. The analyzing device could be used to control some 
prisms which would change the direction and displacement of the beam. The 
analyzing device could be a simple null device if the beam passed through the 
prisms shortly before reaching the analyzer. L. U. Kibler of Bell Telephone 
Laboratories described this idea in an unpublished memorandum in March, 1962. 



2120 THE BELL SYSTEM TECHNICAL JOURNAL, NOVEMBER 1965 




pig i — A guide on a rough surface that produces an irregular distribution 
of forces with a characteristic separation X p . This results in waves with a charac- 
teristic length X„ . 

coordinate axis, x, to the axis of the guide will be called ij(x) and the 
distance from the axis of the guide to the ray or any point in the guide 
will be called 5(x). See Fig. 2. 

For the case of the continuous lens-like medium, suppose that the 
refractive index can be described as a function only of distance, 8, from 
the axis of the guide. A medium free of spherical aberration will be 
defined as one obeying the relation 



n(B) = rc exp (- C8 2 /2) « w (l - C8 2 /2). 



(1) 



C is the specific convergence 3 of the guide, which is approximately the 
convergence, in diopters or reciprocal meters of a one-meter segment of 
the guide if that convergence is small compared to one diopter. 

The equation of the trajectory of a ray in a medium of slowly, smoothly 
varying isotropic refractive index is 



("l) = gH 



— ( n ^r- ) = grad n. 
ds 



(2) 



In this equation ds is an element of length along the path of the ray, 
r is the vector position of a point on the ray and n is the (isotropic) 
refractive index at that point. 



LENS GUIDE 



CONTINUOUS GUIDE 




LIGHT RAY 8(X) 



AXIS OF m 
GUIDE 




Xn = nL 



Fig. 2 — Symbols used to describe a continuous guide (left) and a thin lens 
guide (right). 
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If the slope of the ray, dS/dx, and the slope of the guide axis, dy/dx, 
both remain small, then (1) and (2) may be combined to give the fol- 
lowing differential equation for the trajectory of the ray, 

( ll + ('8 = - ^. (3) 

dx 2 dx 2 ' 

This is like the familiar forced, undamped harmonic oscillator equation. 
For a guide composed of a series of thin, aberration-free lenses of 
equal convergence, c, separated by a distance L, we can obtain the fol- 
lowing analogous difference equation for the trajectory from simple 
geometrical construction. 

cL 

H«n+1 + *n-l) - *» + "2 5 » = -(MZ/n+l + V*-l) - Vn) (4) 

8 n is the displacement of the ray from the axis of the nth lens and y n 
is the displacement of that lens axis from the straight ar-coordinate axis. 
(See Fig. 2.) Note that the product nL is equal to x. Hence, once we 
define a shape, y(x), for the guide axis, the right-hand side of either (3) 
or (4) can be written down. 

III. SOLUTIONS IN TERMS OP GUIDE SHAPE 

Suppose we pretend that the shape of the guide is periodic with a 
length, A, for the periodicity. The path of the axis of the pipe can be 
represented by the Fourier series: 

•«-5'(£M¥-«)- (5) 

The phase factor ipj is used to avoid writing cosine terms in the series. 
We can avoid a constant term (j = 0) by proper choice of the position 
y = 0. Making the substitution kj = 2irj/A, we may write 

y(x) = £ rOW sin (*j*-«). (6) 

3.1 The Case of a Continuous Lens-Like Medium 
By substituting a trial solution of the form 

5(.r) = I) A (kj) sin (kjx - *,) + B(k t ) sin (VCx - fc) (7) 
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into (3), along with the shape spectrum (6), we obtain the results 



A(kj) = k ' Y{ **\ and +j=<pj. (8) 

If we add the boundary conditions 

6(0) = (r/5/r/.r)x=o = (9) 



we find that 



where 



B(kj) = -A(kj)Hm Vj /mx^ (10) 



COt£y = ~hi cot <PJ ■ (U) 

Note that \/C is the natural angular wave number of oscillations of 
the ray about the axis of the guide if the guide is straight. 
Inserting these results into (7) we finally obtain 

5(.r) = S J _ -A cos <pj ( sin kjx - ~^ sin y/Cx\ 

— sin ^(cos kjX — cos y/Cx) . 

At this point in the derivation we introduce the random, statistical 
character of the problem by asserting that the phase factors, tpj , are 
random. Then we can no longer specify y(x) or 8(x) but we can specify 
a mean square value of each and relate one to the other. Since the mean 
square value of a sinusoidal function is half the square of its absolute 
value, we can rewrite (5) as 

(v'(x)) ^ ? = I Z FV*/). (13) 

Similarly, from (12) we obtain 

<**» - '<»> - 1 g (!0w)' [( si - *> - ^ ^) (14) 

+ (cos kjx — cos -y/Cx) 2 . 
If we make the substitution 



e, =-7--l (15) 



Vc 
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and do quite a lot of algebra and trigonometry, we find that (14) can 
be written in the following form: (Remember that e,- ranges from —1 to 
+ ".) 



4 sh\(e J \/ r Cx/2) _ sin tjy/Cx sin 2\/Cx 



in 2 VCr . 



If the length, A, of the periodicity of the shape of the guide is allowed 
to become very large, we may replace the summation signs in (13) and 
(1(5) by integrations over the variable t, and we may drop the index j. 

If c F 2 («) remains finite as e approaches infinity, then large values 
of e will not cause the integral for (5(.r)) 2 to diverge. (Slightly weaker 
conditions could be applied.) 

For large values of x, the first of the four terms in square brackets in 
(l(i) will dominate if F(e = 0) is nonzero. Thus we obtain the following 
result 

Vt \ ~ l f" v*t \ (1 + e)4 4 sin 2 e^/Cx/2 A 



Wc T7 . 



(17) 
F z (e - 0) for large .r. 



The proper normalization of the function F(e) can be obtained from 
the following relation : 

7 =\j_ x Y\e)ch. (18) 

(Compare (13). It is assumed that the relative spectral distribution 
for F(e) is known, and we may also assume that the mean square value 
of y is known. See Section V for examples.) 

3.2 The Case of a Series of Lenses 

An exactly analogous but slightly more complicated procedure can 
be followed for the case of a series of lenses. We can keep the same 
function, (6), to describe the displacements of the lens centers if we let 
.t be an integral multiple, n, of the lens separation, L. (See Fig. 2.) Equa- 
tion (7) for the ray trajectory may similarly be retained with the same 
expression replacing .r. 

However, when (7) is substituted into (4) rather than (3), we obtain 
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a somewhat different expression for A(kj). The phase factors \p„ are 
the same. 

. ,, N (1 — cos kjL)Y(kj) j , , Q /^ 

A(kj) = and fa = tpj • (8 ) 

— — (1 — cos kjL ) 

We use boundary conditions analogous to (9), 

8( n = 0) = 8(n = 1) = (9') 

and we obtain expressions for B and & : 



B(W--A(WS? (10) 

BID i; j 



and 



, fc cos fcL — cos /c c L + sin /cL cot <py (11 '\ 

COt fy = : r—j • Ul ) 

sm k e L 

In (ll'), k c is the natural angular wave number of oscillations of the 
ray about the axis, corresponding to y/C in the case of a continuous 
lens. In this case it is the first positive real solution to the following 
equation, if such a solution exists. 

COSfe c L = l-y (19) 

The denominator in (8') may be written cos kjL — cos k c L. (When 
cL > 4 the lenses are separated by more than four times their focal 
length and they will not guide the beam.) 

The analog of (12) would be like (12) except that terms from (8'), 
(ll') and (19) replace the terms from (8) and (11). It is not worth- 
while to write it out here. 

Equation 19 has an infinite number of real roots if cL < 4. We have 
chosen to describe the ray trajectory using only the first root, but we 
cannot ignore the rest entirely. If the pipe containing the lenses has an 
appreciable amplitude of waviness near one or more of these higher 
wave-number roots, the ray oscillations will grow because the lenses 
will be displaced as if there were additional amplitude in the funda- 
mental, long wave-length component. These higher roots correspond to 
waves in the pipe with length less than 2L or twice the lens separation. 
In choosing the lowest root for k c we are following Brillouin's example 5 
and restricting values of k to the "first Brillouin zone" of the periodic 
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array of lenses. Remember that kj = 2wj/A where A is the length of 
the periodicity in the shape of the guide, which will ultimately be made 
very large. We can rewrite (6) in the following form in order to limit 
the range of fry from zero to ir/L. 

mil °° r /o_i \ 
vim) =EE YV^r + k n sm ( ^ nL + *w.d 

, v f 2ir(h + 1) ,\ ■ n T r 

+ Y I y fcy 1 Sill (fcyWL + ?*./.«) 

(nL is the a; coordinate of lens number n. See Fig. 2.) Then the analog of 
(13) may be written as 

t A/2L (13 ) 

^ y-i 

When written in terms of aS rather than Y, the analog of (14) looks 
like this: 

VM = I E tfflW ( 1 - «" *^ Y 

2 j=i \cos kjL — cos kjjj 

sin fc,L sin k^c 2 



sin A'y.r — 



{(, 

|_\ sin k e L 

+ I cos kjX — cos k c .i 



(cos kjX 

,L)) 2 ]. 



14') 



■ sin fr r .r , . , 
+ . , - (cos k,L — cos k 
sin fc r L 

It is easy to show yourself that if L approaches zero then k r — > -\/C 
and (14 ) approaches (14). 

Next, we make the substitution 

ej= 'it~ l (15 ' } 

in (14 ) and neglect terms of order (eyfr c L) 2 or smaller in the result. 
The neglect of these terms is justified because the result we are seeking, 
as before, turns out to depend only on extremely small values of e 
when x is large. The result of this and a large amount of manipulation 
is the analog of (16) which follows: 
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W SHti) (I - cos k c L s 



fix) 



1 W £fo) A 

2 & /c c 2 L 2 ~ \1 



2 ^1 k?V \1 + cos k c L 

* \yi (16,) 

- 2 sin Avr) - tan * C L f 1 + cos -^ 2 cos k c xj J J 

1 /l - cos fc c L\ A -f L „ 2 / v (u(ej) v\ 
= 2kJD \l + cos W £ * U;; V e/ <;/ " 

As before, we let A go to infinity and change the sum to an integral. 
= 2 . , 1 /l - cos k c L \ f (I/tclH -, v /*(•) , A , 

Although the integral of w/e has a logarithmic divergence at e = 0, 
the divergent part is an odd function of c and the integrand over a small 
finite range around zero is small. The dominant term in the integration, 
for large values of x, comes from the u{e)/t part. Thus we find that 
for large x 

- cos MA /• ( ' / *««- 1 o2/ , 4 sin 2 (efc e ^/2) 



= 2 . , 1 /l - cos k c L\f ( " kch) l #, 

8 ix) *2kjLA l + cos k cL )L S( ' 



2k?U \1 



~ GO8kcL ^.2Tk c xS\e = 0) 



(I/; 



+ cos k c L, 
where, from (13 ), 

s 2 (e = o) = g r 2 (^ /i + fc.) + r 2 (^ (A + D - ft,) . 

Again, when L — > we obtain the same result as from (17). 

IV. THE SHAPE OF THE AXIS OF THE GUIDE 

We might arbitrarily guess a spectral distribution for the Fourier 
components, Y(k), of the shape of the somewhat irregular pipe. It may 
be somewhat more meaningful, however, to guess a spectrum of forces 
on the pipe and to obtain the shape from that. The latter procedure 
enables us to see the effect of stiffness of the pipe, and of the spectrum 
of applied forces, on its shape. For simplicity we will assume that the 
pipe would be perfectly straight, or at least relatively very straight, in 
the absence of the forces that are to be applied to it. 

Suppose the pipe has a modulus of rigidity a = EI where E is Young's 
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modulus and / is the geometric second moment of a cross section of the 
pipe about a transverse axis. The curvature of the pipe is related to 
the local moment of torque, M (:«■), through the equation 1 

^=«M. (20 ) 

(IX- (J 

The second derivative of the torque moment is the transverse force 
per unit length, p, applied to the pipe, so that 

A( -' -) =Vv). (21) 



<Ia 



.i 



a 



If the net force per unit length on the pipe is represented by the expres- 



sion 



p(x) = Z P(*y) sin (kjx - Vi ), (22) 



;=i 



we obtain the following expression for the Fourier transform of the shape 
of the pipe by comparing (0), (21) and (22): 

Y(k) = «*> . (23) 

krtr 

A convenient analytic expression that can be used to approximate a 
reasonable spectrum of forces is the following, where /'„ , k , m and n 
arc adjustable parameters. 

P(k) = Po((A;/feo)7(l + (fc/fco) m ))*. (24) 

In this case (28) becomes 

Y(k) = (Po/ti<r)((k/h) n -*/(l + (A-Ao)'")) 5 . (23') 

Curves of Y and P as functions of k/fa are shown in Fig. 3 for the spe- 
cific values in = 11 and n = 8. The factors outside the square roots are 
omitted. (If necessary we could use a sum of several such expressions 
with different parameters in each without greatly complicating the 
results, but we shall restrict the analysis to one.) 

We will make use of the following formula several times in the follow- 
ing analysis.' 



frfg g - . /(,'+i) ,v (25) 

■»(— i— ) 



Let K represent the ratio />/A n . 
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5.0 



Fig. 3 — Relative amplitudes of sinusoidal components of forces applied to 
the guide (P), and of the shape of the pipe (>'), as a function of relative wave 
number, k/k ■ k c i/k is the ratio of k c to k a used in the examples on the left side of 
Table I, while k c -</k Q is the ratio used in the examples on the right. 



The mean square value of P(K) is 



? 



2 h 



p\K)dK = ^ 



,rPo 2 



r K n 

'o 1 + 



+ K> 



dK 



(26) 



2m sin ((n + 1)t/w) ' 



The mean square displacement on the guide axis from the straight x 
axis is (cf. (23)) 



7 = lf o r {K)dK = \{^f Q J^ m( iK 



(27) 



2m sin ((n — 7)ir/m) ' 

Hence y and p are related through the expression 

./" (w+ IK 

4/ sin 

\ m 



V 

/Co o" 



sin 



n - 7)ir 
m 



(28) 



We will also define characteristic wavelengths, X p and \ y for p(x) 
and y(x) since they are more tangible than the parameter k . The char- 
acteristic wavelengths will be taken as 2tt divided by the angular wave 
number corresponding to the "center of gravity" of the integrands of 
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(26) and (27). Thus we obtain 

Xp = ?I f rr^r-dK I f -Jf—.KdK 

h Jo 1 + K m / Jo 1 + K m 

(29) 
27rsin ((n + 2)7r/w) 



and similarly 



k'o sin ((n + l)ir/m) 



(n - 6)tt 
n sin 

K = £ r- " • (30) 

fco 8 . n (n - 7)* 

m 

For very narrow force spectra X„ and X p are both equal to 2t/Ico , but 
when the spectra are broad, X„ tends to be larger than X p , since short 
wavelength force components are unable to bend the stiff pipe much. See 
Fig. 1. 



V. NUMERICAL ILLUSTRATION 

The following numerical illustration serves mainly to show how the 
preceding results could be used. The parameters chosen are the ones 
that one is likely to know in a real situation. The actual numbers might, 
vary considerably, but it is not hard to estimate how much changes in 
various parameters would affect the results. What is most important is 
to know the order of magnitude of the effects using reasonable param- 
eters. The following example forms the left half of Table I. 

5.1 Shape of Hie Guide Axis 

Let us suppose we have a guide composed of a series of lenses or a 
lens-like medium perfectly centered in a cylindrical pipe. We will sup- 
pose that curvature of the pipe axis is due largely to strain from ex- 
ternally applied forces. For purposes of illustration, we will suppose 
that the main force on the pipe is due to its own weight. We shall sup- 
pose that it is lying on an irregular bed which supports the weight of 
the pipe at more or less random intervals averaging about X p meters 
apart. See Fig. 1. (The spread in the intervals is determined by m and 
win (24).) 

Suppose the pipe is made of steel whose density is p = 7.8 X 10 3 
kg/meter J and whose Young's modulus of rigidity is E = 2.0 X 10 9 
newtons/meter 2 . 



—1 
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Let the pipe be 10 cm OD and 9 cm ID. Let the distance X p be one 
meter. 

Let us use the smallest integer values of m and n that are consistent 
with the conditions on (25) for all the applications of that equation 
(see (26) to (29)). We find that these values are m = 11 and n = 8. 
(This gives a rather broad spectrum of wavelengths of applied force. 
See Fig. 3. Larger values would give narrower spectra, which would be 
appropriate if the pipe were supported at regular intervals.) 

Using these numbers in (29) we find that A = 3.27 meters" 

Equation (30) then gives X„ = 3. 08 meters. 

The stiffness of a cylindrical pipe is related to its OD and ID through 
the equation 

a = -^tf((OD) 4 - (ID) 4 ) (31) 

64 

which gives a = 3.38 X 10 3 newton meters 2 for our example. 

We shall assume for simplicity that the mean square amplitude of 
force variations, p, is equal to the linear density of the pipe times the 
gravity constant, which is correct to within a small numerical factor 
when the weight of the pipe determines all the forces. We thus obtain 

p = p<7((OD) 2 - (ID) 2 ) tt/4. (32) 

For our pipe we get p = 114 newtons /meter. 

Equation (28) then gives us y = 0.407 X 10" 3 meters. 

We shall need to know P when we compute the beam oscillation 
growth rate. Equation (26) gives us P = 222 newtons/meter. 

5.2 Beam Oscillation Growth Rate 

Let us first suppose we have a continuous lens-like guide with spe- 
cific convergence C and the shape defined in part A of this section. We 
use (23') to evaluate Y at k = \/C = 0.5 meters -1 , which is the_value 
at e = 0. If we let C = 0.25 diopters per meter, we get Y{\/C) = 
0.572 X 10~ 3 meters. The ratio \/C/k is 0.153 in this example. It is 
labeled k c i/k on Fig. 3. 

Inserting the value of Y into (17) gives l(x) = y/x- (0.358 X 10" 3 ) 
meters. In 100 meters, the root mean square amplitude of oscillation is 
3.58 mm and in 10 kilometers it is 3.58 cm. 

Next let us consider a series of thin lenses that are separated by twice 
their individual focal lengths, (L - 2/c), and that give the same 
angular wave number for ray oscillations as in the preceding example, 
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k c = 0.5 meters -1 . From (19) we find that cos k c L = 0, which gives 
L = ir or 3.14 meters and c = 2/ir or 0.637 diopters. 

Now we evaluate the sum for S 2 (e = 0) in (17'), using (23'), as in 
the first example. The first term of the sum is equal to F 2 at 0.5 meters -1 , 
whose square root we already evaluated in the first example. The 
next term is Y 2 evaluated at 1.5 meters -1 , the next at 2.5 meters -1 , etc. 
The terms rapidly diminish in size. Equation (17') gives 5 = -y/x- 
(0.830 X 10 -3 ) meters, which is not much larger than the result for the 
continuous lens-like medium with the same ray oscillation wave num- 
ber. 

The preceding examples illustrate the fact that the results are essen- 
tially the same for a series of lenses or a continuous lens-like medium of 
equal k c when the characteristic wavelength, X„ , of irregularities in the 
pipe axis is longer than the separation of the lenses and when the lens 
separation is not too near to the limiting value of four focal lengths. 

The input parameters and results of the preceding examples are listed 
on the left side of Table I. 

The right side of Table I shows what happens when the characteristic 
distance between bumps on the ground, X p , is increased to 4 meters, 
keeping other parameters the same. There is a catastrophic increase in 
rate of growth of beam oscillations because the root mean square am- 
plitude of waves in the pipe is greatly increased while the value of Y at 
k remains near its maximum. See Fig. 3. The ratio \ZC/k or k c /k in 
these examples is 0.611 and is labeled ka/ko on Fig. 3. 

A short FORTRAN computer program is available upon request to 
anyone who may wish to enlarge on Table I using other values of the 
input parameters. 

VI. CONCLUSIONS 

The results are strongly dependent on some of the input parameters, 
but we hope the examples represent the general magnitude of the factors 
one might have to work with if such a light guide were built. By using 
a suitably spaced periodic set of supports for the guide, we could prob- 
ably make Y- or S 2 very small at \/C or k c . The peak or peaks in the 
function Y 2 (k) or S 2 (k) could lie elsewhere. Thus, the contribution to 
wave growth due to forces on the pipe might be considerably reduced. 

The computations did not consider waves in the pipe axis due to 
tolerance limits that would arise in manufacture of the guide or in linking 
sections, but these waves could obviously be incorporated as an addi- 
tional term in Y 2 (k) if they could be estimated or measured. 
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It seems probable that, the beam would have to be recentered at 
rather frequent intervals along the guide unless the aperture were very 
large or the pipe were extremely stiff and straight. 
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